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The dynamics of entanglement and the phenomenon of entanglement sudden death (ESD) [l| are 
discussed in bipartite systems, measured by Wootters Concurrence. Our calculation shows that ESD 
appears whenever the system is open or closed and is dependent on the initial condition. The relation 
of the evolution of entanglement and energy transfer between the system and its surroundings is 
also studied. 
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I. INTRODUCTION 

The dynamics of entanglement in bipartite systems 
has received great attention since the work of Yu and 
Eberly Q, in which the entanglement between the two 
particles coupled with two independent environments be- 
came completely vanishing in a finite time. This sur- 
prising phenomenon, contrary to intuition based on ex- 
perience about qubit decoherence, intrigues great inter- 
ests [2, 0, [H, El, 0, Si. Different from the original work 
of Yu and Eberly [2], in which two independent parti- 
cles embedded into its own dissipative environment and 
there is no any direct or indirect interaction, the effects 
of interaction between the particles and the couplings to 
the same environment have been discussed extensively in 
Ref. Q and [Bj]. In the authors show that for a spe- 
cial initial state the entanglement disappeared in a finite 
time and then revived after a dark period because of the 
interaction between the particles. Furthermore the au- 
thors show that the entanglement sudden death (ESD) 
is sensitive to the initial condition, as proved in [5]. An 
important character in 0, H[ is that the environment is 
dissipative and the transfer of energy between the sys- 
tem and environment is inevitable. Another important 
situation is the dephasing environment, in which energy 
transfer from the system to the environment does not oc- 
cur. Some works have been devoted to this issue @,0]. In 
0], the authors show that disentanglement is dependent 
on the initial condition and temperature of the enviro- 
ment. 

Although the extensive progress in understanding the 
disentanglement, it is still unclear what reason causes 
this phenomenon and what is the physics behind. In 
this paper we try to give an enlightening discussion by 
examining some examples. In the previous studies, the 
couplings with the environment has been considered as 
the critical point of disentanglement and ESD. However, 
as will find in this paper, ESD can also happen in closed 
systems. In our own point, the phenomena of disentan- 
glement and ESD stem from the dissipative terms in the 
Hamiltonian, independent of the coupling with environ- 



ment, and are sensitive to the initial conditions. More- 
over the revival of entanglement can also be explained as 
the effect of the backaction terms in the Hamiltonian. 

The paper is organized as follows. In Sec. II, we 
first review the dynamics of entanglement in (a) Tavis- 
Cumming model [9[ and (b) dephasing system respec- 
tively. In Sec. Ill, the dynamics of entanglement in 
closed bipartite system has been discussed, and ESD and 
revivals of entanglement have also been observed in this 
system. The further discussions and conclusions is pre- 
sented in the final section. 



II. DISENTANGLEMENT IN OPEN SYSTEMS 

The dynamics of entanglement in open systems has 
been discussed extensively p], Q, 0, [H, S 0] • However the 
physics behind the evolution of entanglement is rarely 
touched. In this section we try to present an physical in- 
terpretation for concurrence by examining two examples. 
The first is the Tavis-Cumming model [9], in which two 
non-interaction qubits couple with the same quantized 
field under the rotating-wave approximation. Another is 
the dephasing model, in which two qubits embed into a 
multimode quantized field and the interaction between 
the two qubits is also considered. 



A. Tavis-Cumming Model 

For two spin- 1/2 particles coupled with a single- mode 
cavity field, the Hamiltonian is written as 
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where is bosonic operator and of are the rising and 
lowing operator for spin- 1/2. We focus our discussion on 
the resonating case (c^o = cj), in that the spectrum of 
Hamiltonian can be obtained exactly. For the whole sys- 
tem (system + field), the evolution of the whole system is 
characterized by the interacion between the system and 
field. However from the point of the system (two spin- 1/2 
paricles) , the energy transfer between the system and the 
field happens, which is described by the relaxation term 



aV~ and the backaction term aa^~ . This example is very 
different from Refs [2|, Lfl] , in which the two qubits couples 
with two independent environments separately. In fact 
the couplings with the same field could induce an effec- 
tive interaction between two qubits, which manipulates 
the entanglement between the two qubites Further- 
more, as will be shown in the following, the phenomenon 
of ESD is sensitive to the initial conditions, and the en- 
ergy transfer because of the interaction is directly related 
to the disentanglement and revival of the entanglement. 

The evolution operator U(t) = exp(-iHt) can be cal- 
culated exactly [lOj]. Choose the initial state p(0) = 
(V + r|y?)(y>|)® |0)<0|, in which |p) = sin#|ee) + 
cos 0\gg) with \e(g)) denoting two eigenstates of spin- 1/2 
particles, and |0) is the vacuum state for the quantized 
field. The entanglement in the initial state is measured 
by the purity r and mixing 0. Then the density matrix for 
the system g(t) = Tr E [U (t) p(0)W (t)} can be expressed 
in the basis {|1) = |ee), |2) = \gg), |3) = \eg), |4> = \ge)}, 
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The concurrence c = max{0, £34 — ^£11/022} has 
been calculated, as shown in FigtU for different initial 
states. Obviously the concurrence is fluctuating with the 
rescaled time gt and can be zero in a finite time. We 
have highlight the points with red color that the con- 
currence is vanishing, as the so-called entanglement sud- 
den death (ESD), in Fig IT] in this plot. An important 
point is that ESD is sensitive to the initial state as dis- 
played in Fig HI For the mixed initial state r < 1, the 
ESD happens readily (FigQJ)). However, for the pure ini- 
tial states, the region with vanishing concurrence is com- 
pressed greatly (FigQJi). This phenomenon shows that 
ESD is completely initial-state sensitivity in this model. 

Another important character revealed by FigQ] is the 
fluctuation/revival of entanglement. This phenomenon 
have been discussed qualitatively in Ref. [4]. However 
a clearer interpretation is absent; How the entanglement 
is constructed by the interaction and why and where the 
entanglement decreases with the evolution? What is the 
physics behind the dynamics of entanglement? In this 
part we try to give an enlightening discussion for this 
problem based on this model. Let us first review the 
Hamiltonian Eq. (pQ); Obviously the dynamics of the two 
independent spin- 1/2 particles is determined by the in- 
teraction terms in Eq. (pQ), which is the charge of the 
energy transfer between the system and field. Hence it 
is natural to check simultaneously the variation of the 




FIG. 1: The concurrence for g(t) versus the coupling gt and 
the mixing 0. (a) corresponds to r = 1 and (b) for r = 0.5. In 
both figures the points with vanishing concurrence have been 
highlighted by red lines. 
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FIG. 2: The concurrence (solid line) and Ho (dashed line) 
versus the rescaled time gt. For this plot we have chosen 
(a)r=l, 0=0.2; (b)r=l, 0=1 A; (c)r=0.5, 0=tv/4; (d)r=0.5, 
0=tv/8. One should note that we also plot the concurrence 
below and chose ujq = 1 for plotting Hq. 
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energy (H ) e = (f (*% 
rence in the system. The analytical relation between c 
and (Hq)q is complicated and it is convenient to plot for 
displaying their relation. One could find the variations 
of (Ho) 6 and concurrence is almost in- step, as shown in 
Fig|2l At most cases, with the increment of {Hq) 6 , the 
concurrence increases until the same value with the initial 
states. An exceptional situation appears in Figj2][b), in 
which the second wave peak accompanies with the oppo- 
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site transfer of the energy to the other situations. How- 
ever we could suggest that the dynamics of entanglement 
and its revival in this model stem from the energy trans- 
fer between the system and environment. Moreover we 
also note that the extreme points of the concurrence and 
those of (Hq)q are one-to-one. These phenomena show 
the intimate relation between the concurrence and the 
energy (Ho) 6 in this model. 

However the relation between ESD and the transfer of 
energy is not direct, as shown in Fig 12 Originally ESD 
comes from the cutoff in the definition of concurrence. 
Although the physical explanation for the dynamics of 
concurrence can be available, ESD cannot be explained 
properly based on the transfer of energy in this model. 
A further discussion is needed beyond this model. 



B. Dephasing Model 

In the previous subsection we discuss the dynamics of 
entanglement and ESD in a dissipative two spin- 1/2 sys- 
tem. A novel character is the direct relation between the 
evolution of the concurrence and energy transfer between 
the system and field. Another important situation for the 
open systems is the dephasing case, in which there is no 
energy transfer between the system and environment. A 
typical Hamiltonian for this case can be written as 

H = Hs + He + Hj, (3) 
H E = ^Wjbtbj, 

3 

3 

where two interacting qubits is coupling with a quantized 
environment by Hi . The evolution operator can be given 
by, 

U(t) = exp(-iiJt) 

= e- it ^* + ^i$ b ' ) \ee){ee\ + 

+e -it(n-E^r?M+E i w i o+ + a i+ )| + ^ + | 

in which \g(e)) is the eigenstate of spin- 1/2 particle and 
|±) = (| eg) ± \ge))/y/2 and dj± = bj ± Yj/ujj. Choose 
the initial state p(0) = + r|0)(0|) (g) |0)(0|, in which 
\(j)) = sin 9 1 eg) + cos#|ge), and at time t the reduced 
density matrix g(t) = Tr E [U(t)p(0)W (t)] can be written 
in the basis {|1) = |ee), \2) = \gg), |3) = |+), |4) = |-)} 
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FIG. 3: The concurrence for g(t) versus the rescaled time cot 
and the coupling T/co. (a) corresponds to r = 1 and (b) for 
r = 0.5. In both figures we choose = tt/20 and Q/co = 
3. In (b) the red line high- lightens the points of vanishing 
concurrence. 
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Obviously the off-diagonal terms is damping because of 
the decoherence factor ]\. exp[4fi|(coso;jt — 1)/^|] 

The concurrence for g(t) can be given exactly, as shown 
in Fig J3] with a single mode for simplicity. One notes that 
the entanglement does not decrease, but have a damping 
fluctuation until the same entanglement contained in the 
initial state, as shown in Fig. ©(a). Another facet is 
that when the concurrence is zero for the initial state, it 
is possible for system to generate entanglement, the en- 
tanglement is also a damping-oscillating function of time, 
as shown in Figj3](b). The ESD can happen at some spe- 
cial times. However, for the entangled initial states, there 
is no ESD at any time [l3| ]. 

III. ENTANGLEMENT IN CLOSED SYSTEMS 

In the previous section, we discuss two models coupled 
with a field or environment. In both models the dynam- 
ics of entanglement is dominated by the coupling with 
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FIG. 4: The concurrence versus the time out and (a) the purity 
r, (b) the rescaled coupling J = g/2uo. We have chosen = 
7r/4 and (a) J = 1, (b)r =1/2, and the points with vanishing 
concurrence have been highlighting with red color. On the 
right of the figure (a) and (b), the sectional drawings have 
also provided with the same values of parameters respectively. 
The chosen parameters for the solid, dotted and dashed lines 
are (a)r = 0.35, 1, 0.5 and (b) J = 1, 0.5, 2. 



the field or environment. Next we will show that ESD 
can also appear in closed quantum systems, in which the 
evolution of the system is insulated to its surroundings. 
Let consider two spin- 1/2 particles with Ising-type inter- 
action. The Hamiltonian is 
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The time evolution U(t) = exp(-iHt) is written in the 
basis {|1> = \ee), |2> = \gg), |3> = \eg), |4) = \ge)}, 
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where A 

It is convenient to choose the initial state p = + 
r\(p)(cp\ with \ip) = sin#|ee) + co$6\gg). The entangle- 
ment for p(t) = U(t)pW(t) can be given easily, and the 
concurrence have been plotted in FigJH It is clear that 
ESD happens for some special cases and then the entan- 
glement revives at a later time. We also choose the initial 
state p = + r\(j))((j)\ with \<j>) = sm6\eg) + cos0\ge). 
However there is no ESD . These phenomenon show again 
that ESD is sensitive to the initial conditions. 




FIG. 5: The concurrence (solid line) and Hi (dahsed line) 
versus the time cot. For this plot we have chosen = tv/4 and 
J = 1. The figure (a) corresponds r = 0.5 and (b) for r = 1. 



It is of great interest to check the relation between 
the concurrence and the energy transfer described by 
(Hi)p = (go'A (7 B)p' ft i s easily to note that the con- 
currenc c of p(t) is c ex (Hj) p /J and a figure has been 
drawn for showing their relation , displayed in FigEJ It 
is very interesting to note that the evolution of concur- 
rence and {Hj) p is in- step and there is a one-to-one rela- 
tion of the extreme points between the concurrence and 
(Hi) p . This phenomenon is similar to the case discussed 
in Sec.IIA (see Fig|2]) and show that the dynamics of en- 
tanglement can be related directly to the energy transfer 
in this model. 



IV. DISCUSSIONS AND CONCLUSIONS 

Some conclusions and further discussions should be 
presented in this section. In this paper we discussed the 
dynamics of entanglement and ESD in bipartite systems 
. We found two very different types of the dynamics of 
entanglement. The first case is that there is the effect of 
the dissipation and vice versa, just as shown in Eq. (pQ) 
and ©. Our calculation shows the intimate relation be- 
tween the concurrence and the energy transfer described 
by Ho in Eq. (pQ) and Hj in Eq. (|6|); the extreme points 
of concurrence and the energy of Hq or Hi under the 
time evolution are one-to-one. In the case without ESD, 
the minimums of the energy corresponds to that of the 
concurrence , as shown in Fig. ©(a) and Fig. ©(b). 
However, when ESD happens, the correspondence is de- 
stroyed. 

In fact mathematically ESD stems from the cutoff in 
the definition of concurrence, and so it is difficult to find a 
physical interpretation. We suggest that the concurrence 
could be expressed as c = max{0, /(£?)}, in which f(E) 
is a function of the energy transfer, e.g. f(E) ex (Hi) p /J 
in Sec. III. Then we could define a critical energy E C1 be- 
low which the system must be non-entangled , and above 
which the system must be entangled. When E c is the 
minimal of the energy of the system, then we immedi- 
ately conclude there is no ESD. Contrarily ESD could 
happen. Unfortunately the critical energy may be com- 
plicated and be dependent on the state as shown in Figj2] 
and [5] and prevent us from giving a further analysis. 

Another interesting case is the dephasing model, which 
has been discussed extensively in the past as one of the 
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standard decoherence models [l4| • The evolution of con- 
currence in the model Eq. ([3]) is damping-oscillating with 
time until the same entanglement contained in the initial 
state (see Fig. [3]). Moreover the oscillation of concur- 
rence seems not directly related to the energy transfer 
since there is no energy transfer between system and en- 
vironment. From this phenomenon it may be frustrating 
that the connection of concurrence and the energy is not 
universal. What are the reasons of the dynamics of con- 
currence in the dephasing case is still an open question. 
The phenomenon of ESD does not appear in the dephas- 
ing model when the initial state is entangled. It seems to 
imply that ESD intrinsically originates from the energy 
transfer. However the correct understanding depends on 
the physical meanings of the concurrence. 



In conclusion we have discussed the dynamics of the 
entanglement and ESD in some special models. Some 
interesting phenomena have been presented with enlight- 
ening discussions. In our own points the main obstacle 
of explaining these phenomena is that it is unclear what 
the physical meanings of the concurrence are. Recently 
a great deal of works are denoted to the understanding 
from the energy of the system [15[. However, as shown in 
our paper, it maybe fail when there is ESD happening. 
More recently a physical interpretation of concurrence 
for the bipartite systems has been provided based on the 
Casimir operator in [16[. It maybe opens another door 
to understand concurrence as a physical quantity. 

This work was supported by NSF of China under 
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